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We study the spatial Fourier transform of the spin correlation function G q {t) 
in paramagnetic quantum crystals by direct simulation of a Id lattice of atoms 
interacting via a nearest-neighbor Heisenberg exchange Hamiltonian. Since 
it is not practical to diagonalize the s = 1/2 exchange Hamiltonian for a 
lattice which is of sufficient size to study long-wavelength (hydrodynamic) 
fluctuations, we instead study the s — > oo limit and treat each spin as a 
vector with a classical equation of motion. The simulations give a detailed 
picture of the correlation function G g (t) and its time derivatives. At high 
polarization, there seems to be a hierarchy of frequency scales: the local ex- 
change frequency, a wavelength-independent relaxation rate 1 jr that vanishes 
at large polarization P — > 1, and a wavelength-dependent spin-wave frequency 
oc q 2 . This suggests a form for the correlation function which modifies the 
spin diffusion coefficients obtained in a moments calculation by Cowan and 
Mullin, who used a standard Gaussian ansatz for the second derivative of 
the correlation function. 

The Heisenberg spin chain continues to be of the focus of much experi- 
mental and theoretical study (See RefUand references therein). A number 
of new materials has been synthesized which appear to be nearly perfect ex- 
perimental realizations of the spin-1/2 chain. These experiments have also 
pointed to the need for detailed theroeticals studies of the spin dynamics 
of the Heisenberg chain. In its simplest form the model is defined by the 
Hamiltonian 



where B is the external magnetic field, J is the nearest neighbor exchange 
integral, the second sum is over nearest neighbors. Its low temperature (T <C 
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J) behavior is controlled by the T = quantum critical point and is well- 
described by quantum field theory. The finite temperature dynamics are less 
understood. The infinite temperature regime has been studied theoretically, 
and numerically, but almost exclusively in the unpolarized B = case. In 
this paper, we focus on highly polaralized paramagnet at infinite temperature 
J/T — > 0, but with B/T finite. The 3-d version of the spin- 1/2 system has 
been proposed as a model of spin transport in the paramagnetic hep phase 
of solid 3 HePHl 

Previously, we studied spin transporlPin a 2-d lattice model that demon- 
strated Leggett-Rice type behavior, i.e., characterized by an exchange mean- 
field supporting damped spin-waves. The Leggett equatiorP was originally 
derived for paramagnetic Fermi fluids, but it has been shown that it ap- 
plies to the spin dynamics of any polarized system interacting via quantum 
exchanged Whereas in the liquid longitudinal spin transport is due to par- 
ticle motion and damping of the spin current is due to particle collisions, 
in a lattice of localized spins transport is due to exchange and damping is 
due to spin fluctuations. In the current paper we study 1-d chains and find 
another analogy with the fluid: a relaxation of the spin-current at high po- 
larization, where the spin system begins to acts like a dilute gas with a mean 
free time r. We restrict our attention to longitudinal spin diffusion, although 
a similar description seems to apply for transverse spin transport (the re- 
sults are of a more preliminary nature, however). Specifically, we study the 
polarization dependence of spatial Fourier transform of the spin correlation 
function G q {t) in the hydrodynamic limit q —* 0. The time derivatives of 
this function (or the moments of its Laplace transform) give the diffusion 
coefficient D», and the relaxation time r. At high polarization, there seems 
to be a hierarchy of frequency scales: the local exchange frequency J/h, a 
wavelength-independent relaxation rate 1/t(P) that vanishes at large polar- 
ization P — ► 1, and a diffusive decay rate 7 = D^q 2 . 

Since a lattice of N spins involves 2 N quantum states, it is not practical 
to diagonalize the S = | exchange Hamiltonian for a lattice which is of suf- 
ficient size to study long-wavelength (hydrodynamic) fluctuations. Instead, 
We take a commonly used approach EHSD an d study the S — > 00 limit and 
treat each spin as a vector S(r) with a classical equation of motion, 

|s(r) = 0]TS(r)xS(r') (2) 

r' 

where the sum is over nearest neighbors and £1 = J/h is the exchange fre- 
quency. For convenience, we have moved to the Larmor frame, and have set 
the vector magnitude of the spins |S| equal to 1. Also, we found only quanti- 
tative differences between 1 and 2-d systems, so we have concentrated on 1-d 
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Fig. 1. (a) Spin correlation function for P = 53.7%. (b) Log-Log plot of the 
spin diffusion coefficient vs. 1 — P 2 . form the longtime diffusive behavior 
and the microscopic relaxation time of G'l(t). The last four points have a 
slope of 2.3 ±0.1 



systems, which are computationally less expensive, and for which theoreti- 
cal results are easier to obtain. Although it is knowrPl that the unpolarized 
chain exhibits anomolous diffusion j q oc q a with a slightly less than 2, this 
behavior does not seem to be a present at higher polarizations, where the 
fluctuations essentially pass through one another. 

Longitudinal diffusion was studied by simulating the motion of a hydro- 
dynamic (q —* 0) fluctuation S q (t). We define the correlation function^ 

q{ ' (Sf(O)Sf (0)) 

which is assumed to have an exponential decay oc e- u W q 1 as t — > oo. The 
brackets indicate an average over initial conditions sampled from an equilib- 
rium ensemble. 

We briefly outline the procedure for calculating G q {t) and its derivatives. 
First, an initial spin configuration was obtained by generating a lattice of 
randomly oriented spins Sj with a Boltzman distribution e^ B ' T ' S i . The 
exchange energy was neglected in the Boltzman factor, so the spins were 
initially uncorrelated. The spin configuration was then allowed to evolve 
according to EqEl and the spin correlation function G q (t) was calculated 
from the definition EqEl The average was then performed over N initial 
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configurations where N = 10 6 was typical. Periodic boundary conditions 
were used and the length L of the lattice was chosen to be much greater 
than the "mean free path" given by cr where c = 2£la is the maximum spin 
wave phase velocity, (a is the lattice spacing) and r is the "mean free time" 
(defined below). The equation of motion was integrated with a fourth-order 
Runge-Kutta method, with time step At = 0.1. The typical motion of a 
spin in a time step was 5 degrees or so. The conserved quantities such as 
the total energy, the polarization, and the individual spin magnitudes were 
found to be constant to a high degree of accuracy. 

In simulations to study the short time behavior of G q (t) and its deriva- 
tives the system was allowed to relax for a time greater than the microscopic 
relaxation time r to remove any artifical effects of the uncorrelated initial 
conditions. The higher derivatives were, of course, a good deal noisier than 
G q (t) data and up to N = 10 8 simulations were even averaged in some cases, 
although these were short simulations to study short-time behavior. 
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In Fig. la the correlation function G z q {t) is shown for P = 53.7% and is 
indistinguishable from exp(— 7 9 |t|). From this function we extract the spin 
diffusion coefficient JDm = "y q /q 2 - The polarization dependence is shown in 
Fig. lb, from which we find D|| oc (1 - p2)-2.3±o.l We algo get a detailed 
picture of the short-time behavior of the correlation function. In Figs. 2a- 
d we show the first four derivatives of G q (t). The surprising result is the 
fact that at high polarization the second derivative is not Gaussian but has 
the form of a hydrodynamic correlation function G q (t) = G q (0) exp(— |t|/r), 
where r is a polarization dependent relaxation time that diverges as P — > 1 
and G£(0) oc q 2 . (It is easy to show from Eqs.(2-3) that G£(0) -> -2Q 2 a 2 q 2 
as P — ► 1.) Accordingly, the first derivative relaxes exponentially to its 
quasi-steady value of — "fqG q (t c \), where t c \ is a cut off time r <C t c \ "C 7^ 1 )- 
Thus, we have the relation 

G' g (i cl ) = -D lig 2 G q (t cl ) = fj 1 G''{t)dt = G"(0)t (4) 

giving Dii = G"(0)T/q 2 which relates the hydrodynamic diffusive behavior 
to the microscopic behavior. 

The fourth derivative evolves on an even shorter (in fact the shortest 
possible) time scale ~ I/O. For large polarizations G l q {t) has the form 
form G l q v (t) = G q v (0)f(P,t), where G q v (0) oc q 2 is polarization dependent 
and — ► as P — > 1. The function f(P,t)— > (more precisely to a value 
~ G q (0)/r 2 <C G % q (0)) on a time scale t C 2, where here the cutoff time satisfies 
fl^ 1 < t C 2 < r (See Fig. 2b). Assuming this form for G q v (t) we have the 
relation 

G'£(t c2 ) = q -y^ = j Q Gl v (t)dt = G\ v (0)F(P) (5) 

Putting this into the expression for Du we get our main result 

11 q 2 G q v (0)F(P) { } 

Next we compare our numerical results to the moments method calcu- 
lations of Cowan and Mullin et a/P~EE21 who considered a nearest-neighbor 
pair exchange Hamiltonian to describe the quantum spin-1/2 dynamics of 
paramagnetic 3 He in its hep phase, with a nearest-neighbor pair exchange 
Hamiltonian. Using a Gaussian ansatz for the second derivative they found, 
for arbitrary polarization, ZJy oc [Gg(0)] 3 / 2 /[G* 1 '(0)] 1//2 , which diverges like 
(1 - P 2 )-V2 as P 1. If we use Eq.© and their result that Gf(0) oc 
(1 - P 2 ) as P -»• 1, the diver gence of JDm would be modified, diverging in- 
stead at least like (1 — P 2 )~ 1 . This is still much slower than divergence seen 
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in our numerical simulations. However, the divergence will be faster if the 
integral F(P) — > as P — ► 1, which from Fig. 2d seems entirely possible. In 
any event, our results shift the focus to the calculation of r which depends 
on the short-time (O,^ 1 ) behavior of the fourth derivative. 
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